Laser tunneling from aligned molecules 
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We study multi-photon ionization from N2, O2 and benzene using circularly polarized light. By 
examining molecular frame photo-electron angular distributions, we illustrate how multi-photon 
ionization acts a momentum-selective probe of the local electron density in the Dyson orbitals for 
these molecules. We find good agreement with calculations based on a tunneling model and including 
saturation effects. 



Ionization of atoms, molecules, and solids by strong 
infrared laser fields is a common approach to control and 
study ultrafast phenomena [1, 2^ . Light-matter interac- 
tion in this limit is often modelled as tunneling. For the 
rare gas atoms tunneling accurately predicts the ioniza- 
tion rate over many orders of magnitude [3 . Tunneling 
models place additional limits on the continuum wave- 
function of the photo-electron [4 , but only recently have 
experiments begun to test these more subtle predictions 
[SHi]. Analytical models of atomic ionization suggest 
that the momentum distribution of the bound electron 
is imprinted onto the continuum wavefunction [8j. How- 
ever experiment, simulations, and quantitative theory all 
show the lateral momentum distribution exceeds the an- 
alytical width by 10 - 15% [5]|9l[T0]. Nonetheless, both 
experiment and calculations illustrate that multi- photon 
ionization transfers information from the bound orbital 
to the continuum. 

Work with small molecules is more qualitative. Ex- 
periments show that orbital information, especially the 
location of nodal planes, is transferred to the ioniz- 
ing wavepacket pTHTS] . The wavepacket structure that 
emerges is important for many applications like high- 
harmonic spectroscopy [T4HT6] , laser- induced electron 
diffraction and holography [T7H20], and orbital tomogra- 
phy [2TII22]. Despite its importance for ultrafast science, 
there are few quantitative tests on strong field ionization 
of molecules. In this paper we experimentally study three 
different molecules - N2, O2 and CqHq - to test how well 
theoretical concepts from tunneling of atoms transfer to 
molecular strong field ionization. 

We will show that, for the three molecules we study, 
the saturation intensity depends on the molecular align- 
ment. In the cases of O2 and CqHq the difference is con- 
sistent with the intuitive predictions made by analytical 
theory [8 , however, in N2 it is not. The deviation in N2 
is consistent with prior experimental work [23^ and sug- 
gests modifications to the filtered orbital concept [H [11] . 
We also find the lateral momentum distribution from the 
three molecules is consistent with the values measured in 
atoms at similar laser parameters [5]. 

Our experimental approach is influenced by work on 
atoms. We use circularly polarized light at 800 nm to 



FIG. 1. Concept of the experiment using circularly polarized 
pump & probe laser pulses. The full revival of N2 is sketched 
as the green curve. Just prior to the revival, the molecu- 
lar alignment distribution, sketched in blue, is parallel to the 
propagation axis (k-alignment; inset upper left). A circularly 
polarized probe pulse shown in red singly ionizes the k-aligned 
molecules. For a slightly longer pump-probe delay the molec- 
ular alignment is confined to the plane of polarization (inset 
lower right). We record photo-electron spectra from single 
ionization under both alignment conditions. See the text for 
a discussion of the alignment revivals in O2 and benzene. 



ensure re-collision is essentially eliminated. Thus, the 
momentum along the laser propagation axis (z) is the 
quantum uncertainty of the electron when it departs the 
tunneling barrier. In addition, the photo-electron mo- 
mentum in the plane of polarization allows us to deter- 
mine the laser field strength at ionization [24] . 

We use field-free non-adiabatic alignment to control 
the molecular alignment in the lab frame [25 . This makes 
our approach applicable to neutral molecules without a 
permanent dipole but it imposes a cost of imperfect con- 
trol over the molecular frame. To handle this limitation 
we adopt the method of normalized differences [11] . This 
has the added advantage of removing the need to cali- 
brate the detector response. Additionally, we have used 
a circularly polarized pump pulse to align the molecules. 
In a circularly polarized aligning pulse, the most polariz- 



able axis of the molecule receives a torque into the plane 
of polarization {xy - see Fig.[T]). This creates a rotational 
wavepacket that re-aligns to the plane of polarization ev- 
ery full revival period Tr. For a linear molecule, just 
before the full revival the internuclear axis first aligns 
along the direction of laser propagation, z (Fig. flj inset 
left). We call this "k-alignment" . The molecules continue 
to rotate and a short time later the internuclear axis is 
aligned with the plane of polarization (Fig. IT] inset right). 
We call this "aplanement" . We recorded photo-electron 
spectra under k-alignment and aplanement for all three 
molecules. 

Experimentally, we use 40 fs pulses at 800 nm and in 
the range 4 x 10^^ - 4 x 10^^ W/cm^. The pulses are 
split into a pump and probe using a beam-splitter with 
variable attenuation control in each arm. The circularly 
polarized pump pulse was stretched to ~ 60 fs for the 
small molecules (N2 and O2) or 250 fs for benzene. The 
pump laser pulse is focused by an //25 parabolic mirror 
onto a pulsed, supersonic gas jet (Even-Lavie valve) op- 
erating at 1 kHz. The jet temperature was estimated to 
be 5 K. By exploding the molecules and measuring the 
momentum distributions from N^+ and 0^+ fragments, 
we measured the degree of alignment (cos^ 0) = 0.7 for 
N2 and O2. For benzene we used the proton momen- 
tum distribution measured with respect to the k-axis to 
observe aplanement and anti- aplanement. 

We use a circularly polarized probe pulse to singly ion- 
ize the aligned or aplaned molecules. The probe was fo- 
cused onto the sample of aligned molecules with the same 
parabolic mirror at //1 2. 5. The timing of the probe 
was precisely controlled by an automated delay stage. 
In N2 we used the half revival (aplanement: 4.07 ps, 
k-alignment: 4.40 ps). For O2 we used the quarter re- 
vival for aplanement (2.90 ps) and three-quarter revival 
for k-alignment (8.68 ps). Lastly, for benzene we used 
the full revival for anti- aplanement (87.2 ps) and aplane- 
ment (88.8 ps). In the benzene case the distributions 
sketched in blue in Fig.[l]correspond to the vector normal 
to the plane containing the carbon and hydrogen nuclei. 
Aplanement of benzene corresponds to aligning the nor- 
mal vector along the k-axis (Fig. IT] left). Anti-aplanement 
of benzene corresponds to the normal vector de-localized 
in the plane of polarization (Fig. IT] right). (Since we used 
the full revival for benzene, the timing of the molecular 
aplanement and anti-aplanement is switched with respect 
to the half revival shown in Fig. ITJ) 

The photo-electron spectra were recorded by a veloc- 
ity map imaging spectrometer (VMI) [26 read out by 
a CCD camera (DVC 1312M) and stored on a desktop 
computer. The spectra were corrected for dark noise 
in the camera by subtracting an image taken with the 
laser pulses blocked. Because of the large inhomoge- 
neous electric field in the spectrometer (^ 1 kV/cm) 
and the small mass of electrons, each measured spec- 
trum is a two-dimensional projection of the 3D molec- 
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FIG. 2. Raw lateral momentum distributions for aplaned and 
k-aligned O2 integrated over the plane of polarization, k- 
aligned O2 shows ionization suppression by about 7%. The 
dip at p^ = is due to detector inhomogeneity. Inset: Raw 
2D photo-electron distributions for aplaned (upper) and k- 
aligned (lower) O2. 



ular frame momentum distribution. Each measurement 
contains ~ 2 x 10^ detected electrons acquired at a rate 
of 1.2 electrons/laser pulse. 

An example of the measured raw photo-electron mo- 
mentum distributions from O2 at aplanement and k- 
alignment are shown in Fig. [2J The raw distributions 
are largely a refiection of the strong limits on the elec- 
tron momentum imposed by the shape of the potential 
barrier. Along the z axis the exponential fall off in ion- 
ization probability is due to the increasing tunnel barrier 
in the lateral direction [5]. Tunneling theory predicts 
the lateral photo-electron distribution takes the form (in 
atomic units) [41 [8j 



\i^{pi_)\ (X / dp|| IV^ol exp 






(1) 



where p± is the electron momentum perpendicular to 
the laser field, ?/^o is the bound momentum distribu- 
tion and Ip{P) is the Stark shifted ionization potential 
which depends on /3, the angle between the laser polar- 
ization vector and the major and minor polarizability 
axes [13l |27l l28] . In the plane of polarization (xy) the 
photo-electron drift momentum p = E(ti)/uj where E{ti) 
is the laser field strength at the time of ionization and uj 
is the laser angular frequency. 

The qualitative features in Fig. [2] are largely the same 
for k-aligned and aplaned molecules. This is a refiection 
of imperfect alignment. However, there is 7% less sig- 
nal from k-aligned O2 compared to aplaned molecules. 
This is due to suppressed ionization [29]. From N2 and 
benzene we also observed suppressed ionization from the 
k-aligned or aplaned molecules respectively. The sup- 
pression was 7% was N2 and 5% from benzene. Since we 
will use normalized differences the spectra are not cor- 
rected for the reduced sensitivity of our micro-channel 
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FIG. 3. Experimental normalized difference spectra for the three molecules. Normalized difference is (a — b)/(a + b) where a is 
the photo-electron spectrum measured at k-alignment (benzene: aplanement) and b at aplanement (benzene: anti-aplanement). 
Inset shows the Dyson orbitals of each molecule. Left: N2 at 4.5 x 10^^ W/cm^. Centre: O2 at 2.7 x 10^^ W/cm^. Right: 
Benzene at 4.0 x 10^^ W/cm^ 



plate detector at p = 0. This reduced sensitivity is re- 
sponsible for the on axis minimum and the asymmetry 
in the raw spectra in Fig. [2] By using normalized differ- 
ences, the detector response is calibrated out. 

While this raw width of the lateral distribution from 
the three molecules is largely consistent with the values 
measured in atoms [5 , we are able to precisely compare 
the relative width of the lateral distribution for k-aligned 
vs. aplaned molecules. In O2 the k-aligned distribution 
is 2.3% wider (1/e half width) than the aplaned distri- 
bution, however, in benzene the aplaned distribution is 
12.3% wider. In N2 the aplaned spectrum was 1.3% wider 
than the k-aligned spectrum. 

To observe more clearly the features of the bound or- 
bital in the tunneled photo-electron distributions, we pro- 
cess the measured 2D spectra by taking normalized dif- 
ference, ND = {a — b) / {a -\- b) where a and b are the k- 
aligned (benzene: aplaned) and aplaned (benzene: anti- 
aplaned) distributions respectively. By subtracting and 
dividing two distributions we can partially remove the 
effect of the large contribution from the tunnel barrier 
itself and observe the small modulations due to different 
orientations of the bound state wave- function ipQ. This 
also removes the inhomogeneous detection sensitivity. 

In Fig. [3] we show examples of the normalized differ- 
ence spectra for N2, O2 and benzene. Each spectrum 
shows evidence of the symmetry of the Dyson orbital. 
The Dyson orbital of N2 has maximum electron den- 
sity dit Pz ~ 0. Additionally, ionization perpendicular 
to the molecular axis (from k-aligned molecules) has a 
lower rate than parallel to the molecular axis (aplaned 
molecules) [23]. As we will show, this creates the bright 
red regions at larger longitudinal momentum {px ~ 1.8). 
As one looks at larger lateral momenta {\pz\ > 0.8) the 
ionization probability drops by 10~^ and the signal is 
lost to experimental noise. This is the effect of the tun- 
nel barrier. 



In the O2 spectra, the k-aligned molecules have min- 
imum electron density along the polarization axis. Ion- 
ization along Pz ^ is therefore suppressed compared 
to molecules aligned in the plane of polarization. As we 
will show later, this also means the k-aligned molecules 
will saturate later in the laser pulse when the laser field 
strength is larger. These two combined effects produce 
the four bright corners on the edges of the O2 spectra 
(Fig. [3] center). In benzene there is a nodal plane in the 
Dyson orbital along the molecular plane. The difference 
between the aplaned and anti-aplaned molecules creates 
the regions of positive signal at large lateral momentum 
{\Pz\ > 0.3) shown in red in Fig. [3] (right column). This 
corresponds to the larger electron density in the Dyson 
orbital of aplaned benzene. In Fig. |3] there is a quali- 
tative transition between molecules with maximum elec- 
tron density in the low lateral momentum region of the 
tunnel (N2) and those with minimum electron density 
(O2, CeHe). 

To test the consistency of our observations with the 
tunneling model, we developed a semi-classical numerical 
simulation of tunnel ionization from a Dyson orbital in 
circularly polarized light. To begin, we calculate an ab- 
initio Dyson orbital and Fourier transform it to make a 
3D momentum- space orbital ^p{p). The momentum space 
orbital is rotated into the desired orientation in the lab 
frame. For the anti-aplanement of linear molecules, the 
molecular axis lies along Pz. To calculate an aplaned 
distribution we integrate over 30 angles in the plane of 
polarization (xy). 

The laser field is given by E{t) = 

E e~^^"^^ /'^ [smut I — cosut ]] where r is the 
FWHM pulse duration and i, j are unit vectors in the 
plane of polarization. At each phase of the laser field 
< ujt < 27r the momentum distribution is further 
rotated with respect to the polarization of the ionizing 
field. This makes two assumptions: (1) '^(p) is fixed 
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FIG. 4. Upper row: Ionization probability as a function of time in a circularly polarized 40 fs pulse for k-aligned (anti-aplaned) 
and aplaned molecules. Lower row: Calculated normalized difference spectra (a — 6)/(a + 6) where a is the photo-electron 
spectrum measured at k-alignment (benzene: aplanement) and h at aplanement (benzene: anti-aplanement). Left: N2 at 
4.5 X 10^^ W/cm^ Centre: O2 at 2.7 x 10^^ W/cm^ Right: Benzene at 4.0 x 10^^ W/cm^ 1% Gaussian noise was added 
to each calculated spectrum prior to processing. 



with respect to the molecular frame, (2) the electron 
wave-function adiabatically adjusts to the rotating 
laser field. To calculate the tunneled distribution at 
a given phase, we use Eq. [l] to calculate the lateral 
distribution. The lateral distribution is then mapped 
onto a torus using the classical drift momentum: 
p = ^e"^^"^^^/'^ [cos cjt^ i + sin cjti j] where ti is the 
time of ionization. Integrating over each phase of the 
circularly polarized pulse creates a 3D distribution for a 
particular alignment of the molecule in the lab frame. 

To study the effect of saturation in the resulting photo- 
electron distributions, we track the remaining population 
in the neutral molecule. To calculate the alignment de- 
pendent tunneling rate, we use a modified version of a 
molecular tunneling theory ^30^. We write the instanta- 
neous ionization rate as 



W {E{t), p) = Wadk {E{t), Ipip)) ■ a{/3) 
X / dp_Ldp|| \ipo{p,l3)\ 



^.p Vi^Sl p) 



E{t) 



where Wadk{E) is the atomic tunnel ionization rate de- 
pending on E and Ip{f3) [2 , a{/3) is a free scalar pa- 
rameter and the third factor includes the effect of the 
bound momentum orbital, ipQ. Eq. [2]is similar in spirit 
to many theories on multi-photon ionization of molecules 
which split the rate into an atomic tunneling part and a 
geometrical part [29-31 . For the tt orbitals in O2 and 
benzene a{l3) = const, and Eq. [2] predicts the angle de- 
pendent ionization rate is suppressed along a nodal axis. 
This is consistent with experiment [121 [12 • 



In the case of N2, however, we require a{0) ~ 4a(7r/2) 
to predict the suppression of ionization perpendicular to 
the molecular axis ^23^ . This is true even when including 
angle dependent Stark shifts. 

The ionization probability from the Dyson orbital is 
given by P{t) = 1 - exp [- /_!^ W{E{t'), I3)dt'] and the 

instantaneous ionization yield is dP/dt. We use the in- 
stantaneous yield to weight the contribution from a par- 
ticular alignment angle as a function of time in the 40 
fs pulse. At each time step we average the yield over all 
alignments present in the molecular distribution to calcu- 
late the net yield. From Eq.[2]it is clear the yield depends 
on the alignment of the bound orbital with respect to the 
polarization axis, in addition to the length of the tunnel 
barrier ^ Ip/E. The yield also depends on the remaining 
population in the Dyson orbital. The result is shown in 
the upper row of Fig. |4J 

To compare our tunneling calculations to the experi- 
mental distributions we apply the algorithm to the Dyson 
orbitals of N2, O2 and benzene. In each case, the orbitals 
were rotated to create either a k-aligned/anti-aplaned or 
aplaned set of molecular alignments. We used the experi- 
mental intensity values to calculate the tunneled electron 
distributions and scaled the total number of electrons in 
each calculated spectrum according to the experimen- 
tal values. We did not include alignment distributions 
out of the plane of polarization because we found these 
only change the contrast in the calculated spectra and 
do not yield new qualitative features. To compare with 
Fig. [3] we then created the normalized difference of the 
k-aligned/anti-aplaned vs. aplaned spectra. The result is 



shown in Fig. [4j 

Figure [4] shows the role of saturation in the calcula- 
tion. For N2 (left) and O2 (centre) the ionization rate 
from k- aligned molecules is lower than for the aplaned 
molecules. This means that more population survives in 
the neutral k-aligned molecules later on in the pulse when 
the laser field is larger. Consequently the photo-electron 
distribution is shifted to slightly larger drift momentum 
p = E{ti)/uj. This is the reason for the four bright red 
features on the corners of the O2 distribution (centre col- 
umn in Fig. Is] and [4|. In our calculations saturation also 
plays a role in creating the bright red regions at large 
longitudinal momentum {\px\ ~ 1.8) in the N2 distribu- 
tion. This is again a consequence of the neutral k-aligned 
molecules surviving later into the pulse when the laser 
field is larger. For benzene, the nodal plane along the 
molecular axis means that the aplaned molecules have 
a lower ionization rate and hence saturate later in the 
pulse as verified in Fig. |4] right. The combination of the 
large electron density at non-zero lateral momentum, and 
shorter tunnel barrier around the maximum of the pulse 
envelope produces the bright red vertical stripes along 
the sides of the benzene differences spectra (Figs. [3J|4|. 

The agreement between Figs. [3] and [4] illustrates how 
tunneling probes the local electron density in a static 
molecular orbital. In k-aligned O2 and aplaned benzene 
the Dyson orbital has a minimum in the electron density 
along the polarization axis. In this case we observe that 
bright positive regions in the differences and normalized 
differences appear at larger longitudinal momentum, and 
for non-zero lateral momentum {\pz\ > 0.3). Because 
of the absence of electron density along the polarization 
axis both the total ionization rate and the photo-electron 
amplitude at p^ = are suppressed. In k-aligned N2, 
the Dyson orbital has maximum electron density at zero 
lateral momentum. We observe ionization from those 
molecules is enhanced along the polarization axis and at 
large longitudinal momentum. 

The transition from a local minimum to maximum am- 
plitude along the tunneling axis is potentially a useful 
technique for detecting motion of molecular electronic 
wavepackets. If several ionic states are populated by a 
pump pulse, the ensuing electronic wavepacket could be 
probed in the molecular frame using tunnelling by a few- 
cycle laser field [32]. Modulations in the low lateral mo- 
mentum region of the bound electron density would then 
be imprinted onto the tunnelled photo-electron distribu- 
tion, much like we have shown here for static orbitals. 

By recording MFPADs under different molecular align- 
ment conditions, we have shown how the selectivity of 
laser tunneling is imprinted onto the photo-electron spec- 
trum. Tunneling acts as a momentum selective filter 
that is sensitive to the local electron density. Techni- 



cal improvements in molecular cooling, alignment and 
orientation will lead to improvements in the resolution 
of the tunneling spectra, much like vibrational isolation 
is a prerequisite for high resolution STM. This could 
be exploited in future studies on molecular electronic 
wavepackets [33l [34] . 
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